Phase Modulated Thermal Conductance of Josephson Weak Links 
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We present a theory for quasiparticle heat transport through superconducting weak links. The 
thermal conductance depends on the phase difference ((f)) of the superconducting leads. Branch 
conversion processes, low-energy Andreev bound states near the contact and the suppression of 
the local density of states near the gap edge are related to phase-sensitive transport processes. 
Theoretical results for the influence of junction transparency, temperature and disorder, on the 
phase modulation of the conductance are reported. For high-transmission weak links, D — > 1, the 
formation of an Andreev bound state at e b = A cos(0/2) leads to suppression of the density of states 
for the continuum excitations that transport heat, and thus, to a reduction in the conductance for 
<j> ~ 7T. For low-transmission (D -C 1) barriers resonant scattering at energies e ~ (1 + D/2) A leads 
to an increase in the thermal conductance as T drops below T c (for phase differences near tf> = tt). 

PACS numbers: 74.25.Fy,74.40.+k,74.45.+c,74.50.+r 



The Josephson effect in superconducting weak links 
is perhaps the best known example of macroscopic phase 
coherence. In addition to the superconducting tunnelling 
current, j s = j c sin</>, Josephson showed that the elec- 
trical current through a tunnel junction includes Ohmic 
terms, jd = (co + ci cos<p)V, where <p is the phase differ- 
ence between the two superconductors, and V is the volt- 
age across the junction. These terms describe dissipative 
quasiparticle tunnelling when the junction is biased by a 
voltage. The term, a\Vcos(j), is attributed to "interfer- 
ence" between Cooper-pair and quasiparticle tunnelling 
QI SHI , and averages to zero in voltage-biased junctions. 

Phase-modulated dissipative currents are characteris- 
tic of any type of superconducting weak link. For exam- 
ple, the thermal current through a temperature-biased 
SIS junction is predicted to be a periodic function of <j) 
0. For a stationary phase difference the thermal con- 
ductance is also stationary. Thus, in contrast to voltage- 
biased junctions, the phase-modulated thermal conduc- 
tance does not average to zero. However, less is known 
about phase-sensitive thermal transport across supercon- 
ducting weak links compared with their current-voltage 
characteristics. Recent investigations of heat transport 
through SIS junctions are based on the tunnel Hamilto- 
nian method while Kulik and Omel'yanchuk [|| cal- 
culated the thermal current for the opposite extreme of 
a perfectly transmitting superconducting constriction (a 
"pinhole" ) . In terms of the transmission coefficient of the 
interface potential barrier between two superconductors, 
the SIS junction corresponds to D <C 1, while the pinhole 
corresponds to perfect transmission, D — ► 1. 

In the following we present a theory for quasiparti- 
cle heat transport through superconducting point con- 
tacts for any junction transparency, and as a function 
of temperature and disorder. The thermal conductance 
is sensitive to spatial inhomogeneities of the order pa- 
rameter, particularly changes in phase, which lead to 
branch conversion between particle- and hole-like exci- 



tations 0, and to the formation of low-energy bound 
states in the vicinity of the point contact. The bound 
state spectrum and transmission probabilities for con- 
tinuum excitations are strongly modified by the junction 
transparency, which leads to large changes in the thermal 
conductance of the junction. 

To study quasiparticle transport through temperature- 
biased superconducting weak links, we use the method 
of nonequilibrium quasiclassical Green functions ||. In 
this formalism the advanced and retarded Green func- 
tions, ® A ' R , describe the local spectrum of excitations 
for the system, while the Keldysh Green functions, <8 K , 
carry the information about the nonequilibrium popula- 
tion of these states. Each propagator, & A ' R ' K (p e; R, t) 
is a 2 x 2 matrix in particle-hole (Nambu) space obey- 
ing transport-like equations for excitations of energy e 
moving along classical trajectories labelled by the Fermi 
momentum, pf. We use the notation of Ref. ||, set 
H = l,ks = 1, and consider spin- independent transport 
in spin-singlet superconductors. 

The basic model for superconducting weak links con- 
sidered here is that of a constriction of diameter and 
length on the nanometer scale, much smaller than the co- 
herence length, £ A = Vf/nA, and the bulk elastic and in- 
elastic mean free paths (see Fig. QJ. The potential barrier 
located at z = is characterized by a transmission (or 
reflection) probability, -D(p/) (or R(p/) = 1 — D(pf)), 
for normal-state quasiparticles with Fermi momentum p / 
incident on the interface. The coupling between the two 
superconductors, Si and S2, can then be described by a 
boundary condition connecting the Green functions for 
the two superconductors at the junction interface [{J. 
The order parameter at the junction interface for the 
superconductor Sj is A e l ^ j , and the temperature at the 
junction interface for superconductor Sj is Tj for j = 1, 2. 
The phase difference is denoted by <j> = fa — <f>\ , and the 
temperature bias is ST = T2 — T\ ■ 

Recently, Eschrig || recast Zaitsev's boundary condi- 
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FIG. 1: Temperature-biased ScS weak link geometry. Quasi- 
particle trajectories are coupled by reflection and transmis- 
sion at z = 0. The shaded boundaries define the region 
(~ £a S> a) where superconductivity and the excitation spec- 
trum are strongly modified for cj>\ ^ (j>2- 



tion into a convenient form by parameterizing the qua- 
siclassical Green functions and transport equations us- 
ing Shelankov's projection operators [T^j an d general- 
ized spectral, j R ' A , and distribution, x K , functions obey- 
ing Riccatti-type transport equations. We use the Ric- 
catti formulation of the quasiclassical equations with the 
boundary condition of Ref. Q to solve for the Keldysh 
Green's functions and heat current for superconduct- 
ing weak links driven out of equilibrium by a temper- 
ature bias. The propagators for trajectories incident, 
gR,A,K = feRAKfa . z > Q,z = 0~), and reflected, 



'+ 

R,A,K 



e *t, A ,* = ^r,a,k^ . z < 0,2 = 0") by the interface 
from the S± side can be used to determine the heat cur- 
rent through the interface, 

h = AN f J A eW/ ( T r(^(p / , e )-0^(p / , e ) 

(1) 

where Nf is the normal-state density of states at Fermi 
surface, v/ is the Fermi velocity and A = ira 2 is the cross- 
sectional area of the constriction. The Fermi-surface 
average includes the direction cosine for projection of 
the group velocity along the direction normal to the 



interface; i.e. (...) 

9 = arccos(vj • z). We consider the case in which 5*1 and 
S2 are identical s-wave superconductors with isotropic 
Fermi surfaces; it is straightforward to generalize the re- 
sults to two different s-wavc superconductors. We also 
consider small junctions, a <C £a! in this limit the Joseph- 
son supercurrent through the contact is small by at least 
a factor a/£ A compared to the bulk critical current, so 
pair-breaking corrections to the order parameter by the 
supercurrent can be neglected Thus, to a good ap- 
proximation the order parameters take their bulk equilib- 
rium values, but with local values for the phase. Also for 
small constrictions the thermal resistance of the junction 
is much larger than in that in the bulk, so the tempera- 



ture drop occurs essentially at the junction. 

With these considerations, <S± can be calculated as 
follows. First, we express the Green function in terms of 
Riccati amplitudes For example, 
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N R N A \ if rf - x? 7l A x? + xft R v A 

(2) 

where NfNf = (1 + if f f ; )(1 + J A T A ). The advanced 
amplitudes arc obtained from retarded functions using 
the symmetry 7^ = — (7^)*. In the point contact geome- 
try, the amplitudes and distribution functions for incom- 
ing quantities (lower case) take their local equilibrium 

iA /{e R + i ^A 2 - (e R ) 2 ) , 



values: 7^ 



xf = if ( 
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tanh(e/27j), for j = 1,2. 



Using the boundary conditions of Rcf . |8j| , we construct 
the corresponding functions for outgoing trajectories (up- 
per case), 



R(l + 7272)71 + D(l + 7172)72 

1 + i?7272 + -D7172 
R\l + 7272| 2 2i + D\l + 7i72| 2 S2 — RD\yi 
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where we omitted the superscripts. Inserting these ex- 
pressions into Eq. J2J), and performing the analogous 
calculation for we obtain explicit expressions for 

the phase-sensitive thermal conductance, defined by I e — 
—K((f),T) ST, in the limit ST -> 0. The general result 
can also be used to calculate the heat current beyond 
the linear response limit. A more detailed derivation and 
discussion will be presented elsewhere ^3 ■ 

In the clean limit, e R — > e + i0 + , the thermal con- 
ductance is expressed in terms of the transmission of 
bulk excitations of energy e > A and group velocity, 
v g (e) = Vfy/e 2 — A 2 /e through the junction, 



T)=A deN{e) [ev g (e)} S(e, </>) 



dT 



(4) 
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1) 1 is the Fermi function, Mie) = 



where /(e) 

N f e/y/e 2 - A 2 is the bulk density of states (DOS), and 
S(e, 4>) = D cc (e, </>) + S h(e, 4>) is the transmission coeffi- 



wherc cient for the heat current, which is the sum of 



£ ee (e,0) = D 



(e 2 -A 2 )(e 2 -A 2 cos 2 §) 



e 2 - A 2 (l-L>sin 2 |) 



(5) 



the transmission coefficient for electron-like (hole-like) 
quasiparticlcs remaining electron-like (hole-like), and 



D(l-D) 



(e 2 - A 2 ) A 2 sin 2 I 



e 2 - A 2 (l-L>sm 2 §) 



(6) 



the transmission coefficient for electron-like (hole- 
like) quasiparticlcs with branch conversion to hole-like 
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FIG. 2: a) The thermal current as a function of (f> and barrier 
transparency, D. The thermal conductance is normalized for 
each D by its value at cj> = 0. b) Non-monotonic oscillations 
of K((f>) for D ~ 0.5. c) Normalized transmission coefficient 
at cj> — n showing the resonant transmission for C < 1, and 
the suppressed transmission for D — * 1. The resonance peak 
is at e ros ~ A(l + D/2). 



(electron-like) quasiparticlcs. Here we neglect the angu- 
lar dependence of the barrier transmission and reflection 
probabilities [l5j . 

The heat current density carried by bulk quasiparti- 
cles of energy e in the superconducting leads reduces to 
the normal-state current density, J\f(e)[ev g (e)] = NfCVf] 
the increase in bulk DOS is compensated by the reduc- 
tion in the group velocity. For = 0, D(e, </>), reduces 
to the barrier transmission probability for normal-state 
quasiparticles, D(e, 0) = D. Thus, for = the thermal 
conductance of the ScS contact is simply reduced by the 
opening of the gap in the quasiparticle spectrum. How- 
ever, for (j> the transmission coefficient, D(e,4>), in- 
cludes the modification of the local DOS near the contact 
by the formation of an Andreev bound state (ABS) with 
energy below the gap edge, as well as the particle-hole co- 
herence amplitudes which alter direct (ee) transmission 
and generate branch conversion (eh) scattering. 

The relative importance of the direct (ee) and branch- 
conversion (eh) processes to the phase dependence of 
the thermal conductance depends on the barrier trans- 
parency D. For both processes the ABS plays a central 
role in controlling the phase modulation of the conduc- 
tance. The bound state leads to a reduction in the local 
DOS near the contact and to a corresponding suppression 
of the transmission coefficient for excitations with e = A. 
For moderate to high transmission barriers (D > 0.5) 
this leads to suppression of the thermal conductance for 
<p ?s 7r. For low transmission barriers (D <C I ) multiple 
Andreev reflection leads to a shallow bound state just be- 
low the continuum edge at = A^J 1 — Dsin 2 (|). The 
spectral weight of the ABS is derived from the contin- 



uum states near e = A, which suppresses the divergence 
at e = A at the cost of a large, but finite, resonant en- 
hancement in the transmission of quasiparticles at ener- 
gies, e w A(l + iL>sin 2 (f )), above the gap (Fig. Hi). 
The resonance generates a strong enhancement of the 
thermal conductance as the phase is tuned to <fi = 7T - 
These features, as well as the evolution of the phase- 
modulation of the conductance with barrier transmission, 
are shown in Fig. where we plot the normalized con- 
ductance, k(4>)/k(4> = 0) for < D < I and T = 0.5 A 
(T = 0.72 T c ). For intermediate values of the barrier 
transparency, D ~ 0.5, the phase dependence of the con- 
ductance is a non-monotonic function of <p for < </> < 7r 
(Fig. GJ>), although the amplitude of these oscillations 
is very small. For D = 1, transmission with branch 
conversion drops out (2) h = 0) and the transmission 
coefficient for the quasiparticle heat current reduces to 
D(e,0) = (e 2 - A 2 )/(e 2 - A 2 cos 2 (|)), with a resulting 
thermal conductance in agreement with Ref. [|| for a 
pinhole. 

We compare Eq. 0|for Del with the heat current ob- 
tained by perturbation theory from a tunnelling Hamil- 
tonian (tH) description of SIS tunnel junctions. Based 
on the tH method Guttman et al. @ obtained a heat 
current of the form: Iq + 1\ cos <j). If we expand 5) to 
leading order in the barrier transmission probability, D, 
we obtain the tH result for the conductance from the 
linear response limit of Eq. (2) of Ref. 0b, 



K tH = ANfVf D 



de 



e 2 — A 2 cos ( 
: e 2 - A 2 



dj_ 

dT 



^ • (7) 



This result has an unphysical divergence due to the sin- 
gularity at e = A. In the tH method the divergence is 
regulated by an ad hoc procedure. Guttman et al. re- 
quired a finite temperature difference and A(Ti) ^ A(T2) 
for the gaps of the two superconductors. However, Eqs. 
14101 for arbitrary transparency show that the unphysical 
divergence is not a singularity of the limit ST — > 0, but a 
failure of perturbation theory in the tunnel Hamiltonian; 
which does not include the change in the spectrum near 
the contact. The bound state and the resonance regulate 
the singularity obtained in perturbation theory and lead 
to a thermal conductance that is nonanalytic, ~ DlnZ), 
but vanishes for D — ► 0. The result for k((/)) to order D 
also has non-perturbative corrections to the phase mod- 
ulation of the conductance, which includes terms cx cos <fi 
as well as sin 2 ^ In (sin 2 w), 



K((j)) = K 



• 2 9 i / • 2 <A\ 

Kisin - ln(sm —J 



k 2 sin - , 



(8) 



fcsech 2 (A/2T), k = A N f v f D(A 3 /4T 2 ), 



where K\ 

kq = k dx x 2 sech 2 (xA/2T) J «2 = 
+fc[(4T/A)(l - tanh(A/2T)) 
[1 +^ (A/T)\/x 2 
sech (Ay/x 



and c 



(1 + 1hD)ki 
2 Jg 00 dx xlnx 
Ttanh(A v / ^ 2 TT/2T)] (x 1 + l)" 3 / 2 
1/2T). The relative magnitude of the 
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FIG. 3: The temperature dependence of phase-modulation of 
the thermal conductance for = 7r, normalized by the con- 
ductance at T c , km — ^j-ANfVfT c D. Shown for comparison 
is the conductance for (f> — 0. 



non-perturbative correction, K1/K2, is approximately 
25% for D = 0.01 at T/T c = 0.8 and increases with 
decreasing temperature. 

For general transparency the magnitude of the phase 
modulation of the thermal conductance is a maximum 
for 4> — 7r, except for a small range of barriers with 
D as 1/2. The temperature dependence of the conduc- 
tance for (j> = 7r is plotted in Fig. [3J The thermal con- 
ductance for <j> = 0, is also shown for comparison. For 
moderate to high-transmission junctions (D > 1/2) the 
thermal conductance for <f) ^ is suppressed relative to 
the conductance at cj> = at all temperatures. However, 
for low-transparency junctions (D <C 1) resonant trans- 
mission of quasiparticles just above the gap edge leads 
to increase in the conductance when T drops below T c . 
This effect is pronounced for junctions with D < 0.2; its 
observation would provide a test of this theory of phase- 
induced resonant transmission of quasiparticles. 

The results above were obtained in the clean limit 
for the superconducting leads. In the diffusive limit, 
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<C I <C £oi where I is the elastic mean- free path, 
the excitations and pairing correlations are governed by 
Usadel's diffusion equations for the Fermi-surface av- 



eraged propagators, 



/dW^CP/jR). For 



an ScS contact Nazarov derived a boundary condition for 
the propagator in diffusive conductors |l4| . 

2D [0 2 , fli] 



1 



(9) 



D 



+ D({q 2 ,9i} - 

where Qj is the Keldysh matrix representation for the 
0* ,R ' K , i?b is the barrier resistance for normal leads and 
(. . .} d = J dDp(D) . . . I J dD p(D) D is an average over 
a distribution of channels with transmission coefficient 



D characterizing the interface. For a single channel con- 
tact with transmission coefficient D application of Eq. 
El yields the result from Eqs. 14161 obtained in the ballis- 
tic limit with the replacement: NfVfAD/4 — > l/2e 2 i?b- 
Thus, the phase modulation of the thermal conductance 
of small Josephson weak links are the same in the clean 
and diffusive limits. This result is due to the cancel- 
lation of impurity renormalization of the diagonal and 
off-diagonal self-energies in the propagators for s-wave 
superconductors up to order a/£ A , and that Nazarov's 
boundary condition is based on a junction model with a 
central layer and interface described by Zaitsev's bound- 
ary condition. 

In summary, we have presented a theory for heat trans- 
port through Josephson weak links. For high transmis- 
sion junctions the reduction in states with e > A, result- 
ing from the formation of an ABS near the Fermi level, 
leads to a suppression of the conductance near <f> = -k. 
For small transparency, the presence of a shallow bound 
state produces a resonance in the continuum just above 
the gap edge. This leads to an increase in conductance as 
the temperature drops below T c for junctions with <fi as n. 
For a single channel contact, these results are insensitive 
to impurity scattering and hold in the clean and dirty 
limits. 
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